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1. Introduction . Spacecraft observations have shown that ions are 
routinely accelerated to energies from - 10 keV to - 10 MeV at inter- 
planetary shocks (1). Of particular relevance in understanding how such 
ion events are produced are the angle 0 between the mean shock normal 
and mean upstream magnetic field_and thelevel of magnetic fluctuations 
in the shock's vicinity. When 0^ < 90° (i.e., a nearly perpendicular 
shock) and the magnetic field quasi-laminar , pitch angle scattering is 
infrequent, and ions are accelerated primarily via the scatter-free 
shock drift process by drif ting_along the tj x % electric field parallel 
to the shock surface (2). When 0 g > 0° (i.e., a nearly parallel shock) 
and the level of magnetic wave activity is high, pitch angle scattering 
is frequent, and ions are accelerated primarily via the diffusive shock 
acceleration process by being compressed between converging magnetic 
irregularities fixed in the upstream and downstream flows (3). 

Although a small minority of ion events observed near 1 AU 
have features indicative of the acceleration process being predominantly 
either scatter-free shock drift or diffusive, most events are more 
complex, and are associated with shocks having intermediate values 
°f 0 Bn and various levels of wave activity. To model these cases, we 
have developed a test particle simulation that integrates along ion 
orbits in a system described below. 

2. Model . Let K[X,Y,Zl denote a system fixed with the shock, with the 
shock lying in the Y-Z plane, and the unit vector X = -n (n - shock 
normal), so that X K 0 upstream (subscript 1) and X > 0 downstream 
(subscript 2). The vectors tj^ = U^ (cos 6^, 0, sin 6,) and 1L. = 
Bq^cos 0. , 0, sin 0^_) denote the upstream plasma velocity and mean 
magnetic field (0 { = 0 fin ) , respectively. In K the U x lL electric field 
is along Y. Given in addition the upstream Alfven Mach number and 
plasma beta p j , the mean downstream conditions are obtained by solving 
the MHD jump equations (with a ratio of specific heats of 5/3). 

We assume that the injected ions are true test particles 
(i.e., they are not coupled self -consistently to the shock and/or 
waves), and neglect shock structure by requiring (ion gyroradius) » 
(shock transition). Pitch angle scattering is introduced by superposing 
upon Bq (i - 1 or 2) a zero-mean, random magnetic field component *b^(z) 
which, in either the upstream or downstream plasma frame, varies only 
with coordinate z along 1$^ , is transverse to^ Qi , and is static (i.e., 
scattering is elastic in either plasma frame). The field 1>^(z) is a 
superposition of N circularly polarized plane waves with wavevectors 
along Bq.. The amplitude of each Fourier component with wavenumber k is 
derived from a power spectrum P(k) using a well-known technique (4). 
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The simulation proceeds first by constructing thet>^(z), then 
by solving the Lorentz force equation for the particle orbit using the 
field IS. (z) = iL, + ^(z) in the appropriate plasma frame, performing a 
Lorentz 1 transformation between plasma frames at shock crossings, and 
continuing until pre-set spatial and/or temporal boundaries are crossed, 
whereupon the process is repeated for a new particle* 


3* Results* Acceleration at interplanetary shocks was simulated using 
the form P(k) in Figure 1, which shows sample upstream spectra with Bq = 
B Q1 = 5 nT. Wavenumber k in the solar wind or plasma frame (top and 
right axes) was related to frequency f in the spacecraft frame (bottom 
and left axes) via k - 27tf/V gw for a solar wind speed V g ^ = 4 x 10 
cm/s. Spectrum A (correlation length 2 x 10 1 cm) represents that of 
ambient transverse Alfvenic fluctuations in the interplanetary medium 
(5). Spectra B and C (correlation length 1.3 x 10 8 cm) represent trans- 
verse MHD waves presumably driven by ions streaming upstream from the 
shock (6,7). To model the observed spectra, we damped the upstream 
shock-associated field amplitude with the function [1 + |xj/Xi*] A , 

where Xi* = Di cos 0,, with Di measured along Bq. from the shock. For 
- 1.5 x 10 fl cm = 0.01 AU (» 1.3 x 10 8 cm), tne integrated power or 
variance cr of C at |x| = X^* is 0.06 B , 1/4 that of spectrum B at X = 
0, where a 2 ■ 0.25 B 2 . Downstream spectra used were similar to those 
in Figure 1, except ^ - 1.5 x 10^ cm = 1 AU, and a =0.36 B^ to 
model the enhanced wave levels observed downstream of shocks. Realiza- 
tions of spectra A and B for N = 4096 were superposed to form S.(z), 
with linear interpolation used between grid points during orbit integra- 
tions. The upper scale in Figure 1 shows the resonant proton energy at 
which proton gyroradius p ~ k~*. 


Figures 2 and 3 show re- 
sults for 6 ^ = 0° , 0 ^ = 60°, = 4 

x 10 7 cm/s, B 01 = 5 nT, = 

8, = 1, and protons injected up- 

stream with energy Eq = 10 keV in 
the upstream plasma frame. For this 
case, Bo 2 /Boi - 3.2 and the plasma 
density jump r = 3.7. We define the 
scale time = eBQ^/mQC = 13 sec 
(upstream proton gyroperiod) and 

scale length = V 0 T 0l^ % = x 

10 8 cm (upstream gyroraaius of a 10 

keV proton with speed v Q = 1 .3 x 10 8 
cra/s). 


In Figure 2(a) kinetic en- 
ergy E/Eq in the shock frame is 
shown versus distance X/p Q . from the 
shock (upstream at left, downstream 
at right of X = 0), In 2(b) and 
2(c) E/Eq and X/p Q1 are shown versus 
time t/x^. during UA the orbit. Per- 
iods A-Hin 2(a) are marked atop 
2(b). The largest and most rapid 
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Fig. 1. Wave power spectra 
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energy gains clearly occur during the shock drift phases of acceleration 
(e.g., = 8 to 9, 11 to 14, 24 to 30, and 40 to 42). In compari- 

son, diffusive energy gains due to compression between the inflowing 
upstream plasma and the shock are relatively small (e.g., t/t m = 14 to 
24 and 30 to 40). The particle in Figure 2 crossed the shock 1 42 times 
and was accelerated from 10 to 560 keV in 43 t.. = 9.4 min. This orbit 

shows features typical of the case = 60°. U1 

In Figure 3 we show the energy spectra that result after an 
elapsed time of 300 t~. * 66 min for 2900 injected protons. The 
quantity Af(E)/AE is the fraction of particles with energy E with- 
in AE centered at the logarithmically-spaced plot points. The upstream 
(solid) and downstream (dashed) spectra denote sums over all particles 
with X < 0 and X > 0, respectively. Error bars are two statistical 

standard deviations in length. 

4 . Discussion . The most striking aspect of Figure 3 is that well- 
formed energy spectra spanning more than two decades in energy are 
produced within only one hour after injection and that these spectra 
bear remarkable qualitative and quantitative resemblance to observa- 
tions. For example, the downstream spectrum from 10 keV to 1 MeV is 
well-described by a power law with spectral exponent ~ 2.2, which is 
quite close to some observed values (1). In addition, the upstream 
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Fig. 2. Sample proton orbit 
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spectrum folds over with decreasing 
energy, similar to observed events 
(8). The cutoff time as well as the 
upstream escape boundary at X * 0.5 
AU produce the steepening of the 
upstream spectrum beyond ~ 1 MeV. 

We emphasize that the 
spectra in Figure 3 are primarily 
the products of the shock drif t 
acceleration process, with pitch 
angle scattering simply providing 
the means to return particles to the 
shock for multiple drift accelera- 
tion phases, as in Figure 2. Dif- 
fusive acceleration contributes 
little to acceleration in this 
case. To estimate the maximum 
effect expected from diffusive 

’, and 
fixed 


Energy (keV) 



Fig. 3. Proton energy spectra 


/Bqj = 1 and r 


acceleration, we set 0^ * 0* 

keeping all other parameters 

as in the 0. = 60° case, injected particles at 10 keV (Bq 2 > wx 
= 3.85 in this case). After 300 the spectra extended no further 
than E/Eg = 10 or 100 keV, illustrating, as one might expect, that 300 
gyroperiods is far too short a period for diffusive acceleration alone 
to produce an extended energy spectrum. Examples of energy spectra for 
various values of 0, from 0° to 75° are shown by reference 9. 


5. Conclusions. We have presented an example from a test particle 
simulation designed to study ion acceleration at oblique turbulent 
shocks. For conditions appropriate at interplanetary shocks near 1 AU, 
we have found that a shock with 0 = 60° is capable of producing an 

energy spectrum extending from 10 1ceV to ~ 1 MeV in ~ 1 hour. In this 
case total energy gains result primarily from several separate episodes 
of shock drift acceleration, each of which occurs when particles are 
scattered back to the shock by magnetic fluctuations in the shock 
vicinity. 
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